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I. INTRODUCTION 

Crystal defects are physical entities of primary importance in solid state 
physics. Owing to the breaking of translational symmetry, they originate local- 
ized levels capable of trapping electrons or holes. On this ground, the physical 
properties of materials such as semiconductors can be modified by suitable dop- 
ing determining the sign of majority carriers. Impurity or intrinsic defects may 
originate localized levels allowing recombinations of electron-hole pairs injected 
into bands by some excitation mechanism, thus affecting the lifetime of free car- 
riers. The content of intrinsic defects at (lattice) equilibrium is ruled by general 
laws 1 . But in reality it strongly depends on the history of the sample dealt with. 
Depending on the nature of the defects, three basic schemes are to be consid- 
ered. They are: the Shon Klasens (SK) scheme, in which a conduction band 
electron recombines with a hole kept in a localized level, the opposite Lamb- 
Klick (LK) scheme, in which a valence band hole recombines with a localized 
electron, and the Prener Williams scheme, in which both electron and hole are 
localized within a two level center 2 . 

At equilibrium, the occupation of electron (hole) levels is ruled by the 
Fermi-Dirac distribution f n (E) (f p (E) = 1 - f n (E)), that is, f n (E) = 
1/ {1 + exp \(E — Ef) /kT]} where E and Ef stand for the actual and the Fermi 
levels. The equilibrium densities of conduction band electrons and valence band 
holes are 3 n = N n exp [(E F - E c ) /kT] and p = N p exp [(E v - E F ) /kT], re- 
spectively, where Eq and Ey are the energy levels of the bottom of the conduc- 
tion band and of the top of the valence band, respectively, N n = 2 (2nm e kT / h 2 ) 3 ^ 2 
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and N p = 2 (2irm p kT / h 2 ) , m e and m p the electron and hole effective 
masses, respectively . It is to be pointed out that distinction between electron 
and hole traps is considered for convenience since only the capture cross section 
determines the trapping properties. Often, in the physics of semiconductors 
a hole trap is presented as a deep electron level, while a hole trapping level 
is presented as an acceptor level. In the physics of ionic crystals it is usual 
to define the trapping properties by means of the defect-charge states. Thus, 
negative charge-defects are traps for holes (positive carriers). On the contrary, 
positive-charge defects are traps for electrons (negative carriers) . But there are 
cases where these rules do not hold ( for example, the U~ centres) 4 . Also the 
distinction of carrier capture for trapping or recombination may be considered 
as conventional, since in both cases we are dealing with changes in level occu- 
pancies. Thus, for a general discussion about capture processes it is convenient 
to follow the model of trap classifications used by Simmons and Taylor (ST) 5 : 
"when a trap is empty it is ready to receive an electron, and thus it is operating 
as an electron trap. When the trap contains an electron, it is ready to receive 
a hole, and hence is a hole trap. (We are assuming that the traps are monova- 
lent.) It is convenient to assume that the traps existing below the equilibrium 
Fermi level are neutral when filled with an electron and that the traps positioned 
above the equilibrium Fermi level are neutral when empty. (...) Thus a trap po- 
sitioned above the equilibrium Fermi level is neutral when acting as an electron 
trap and negatively charged when acting as a hole trap. On the other hand, a 
trap positioned below the equilibrium Fermi levels is neutral when acting as a 
hole trap and positively charged when acting as an electron trap. " 

Departure from equilibrium can be obtained by different excitation sources. 
In this paper, we are considering ionizing photons causing band-to-band transi- 
tions. As the excitation source is turned on, carriers injected into the conduction 
and valence bands are drawn by the several processes occurring in the material, 
including trapping in metastable levels. The latter process allows crystals to 
attain excited states which can be held after the excitation source is turned 
off. In this case, the fundamental state can be reached by thermally stimulated 
processes 6 . Under steady excitation, the population of electron- or hole- levels 
attain a steady distribution after a time which, in some cases, may be very 
long 7 . The statistic of occupancy of the traps may be obtained from two dif- 
ferent points of view. The first of these, which was used by Shokley and Read 
(SR) 8 , considers the rate equations for the conduction and valence bands. The 
second, which was used by ST 5 , considers the rate equations for a particular 
trapping centre. In the non-degenerate case, both these approaches lead to an 
occupation f(E) of trap levels with energy E and density N t (E) which is given 
by 

f(F] c n n + CpNp exp [(E v - E) /kT} 

S[ '~ c n {n + N n exp [(E - E c ) /kT}} + c p {p + N p exp [{E v - E) /kT}} 

(1) 

where c„ and c p stand for the electron and hole capture probabilities, respec- 
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tively. The basic assumption of the two approaches is that, at the steady state, 
the population of electron- or hole-levels can be described by means of a Fermi- 
Dirac-like function with suitable quasi- Fermi levels (QFL). For free carriers it 
was found that 

^=exp[(F„-F p )/fcT] (2) 
np 

where F n and F p stand for QFLs of free electrons and holes, respectively. 

QFLs play a role in the rate of recombination processes. A suggestive idea 
introduced by SR considers the recombination rate at the steady state as a cur- 
rent passing through a resistance , called a "recombination resistance" (RR), 
depending on the kinetic parameters of the process dealt with. In the sim- 
plest cases, recombinations follow from two capture processes, that is, hole and 
electron capture, so that two resistances are to be considered , that is, near 
equilibrium 

R n = kT/nptCn (3) 

for electron capture and 

Rp = kT/pntc p (4) 

for hole capture, n t — ftN t and p t = (1 — ft) N t standing for the equilibrium 
densities of trapped electrons and holes, respectively, and f t for the level occu- 
pancy as given by the Fermi-Dirac function. By presenting QFLs as " voltages" , 
SR proved that in the near-equilibrium approximation the above resistance def- 
initions allow for the equation 



v (R n + R P ) = F n -F p (5) 

which has the formal structure of Ohm's law. SR also derived an equation for the 
lifetime of free carriers, that is, r = rvpRjkT (n + p) where R — R n + R p . About 
this point they concluded that " the effect of a number of different sorts of traps 
may be considered on the same basis. For each variety, the recombination is 
represented by a pair of resistances in series and these series pairs are combined 
in parallel for the entire system" . 

At the level of the SR treatment we cannot immediately give a thermody- 
namic meaning to the product vR 2 which, in the theory of electricity, corre- 
sponds to the dissipated electrical power. In all probability, for this reason the 
concept of RR has remained an unexplored minor outcome of the SR statistics 
although sometimes it is recalled in papers devoted to studies on the kinetics of 
electronic processes in semiconductor based devices 9 . The scope of this paper is 
thus to investigate the recombination processes at the steady state from the gen- 
eral point of view of the thermodynamics of non-equilibrium. To this purpose, 
we remain in the range of linear irreversible processes, where it is possible to ap- 
ply successfully the principle of minimum rate of entropy production (MREP). 
This approach resembles the one used (for a different scope) by V. Maxia 10 . 
The goal of this paper is to give a thermodynamic basis to the concept of RR. 



3 



It is shown that affinities are more suitable to represent "voltages " than QFLs. 
What is better, it is shown that the analogue of "dissipated power " is closely 
related to the rate of entropy production of processes dealt with. Thus, the 
paper is structured as follows. In Section II, basic concepts of non-equilibrium 
thermodynamics as well as the MREP principle are recalled. A suitable vari- 
ational calculus based on the MREP principle is presented. In Section III, 
the results thus obtained are applied to typical recombination processes. Final 
remarks are given in Section IV. 



II. THERMODYNAMICS OF THE NEAR-EQUILIBRIUM 

STEADY STATE 



Non-equilibrium thermodynamics lacks a constructive criterion such as the 
one given by maximum entropy for the equilibrium state. The latter provides 
a starting point for the application of statistical mechanics and thermody- 
namics. However, non-equilibrium thermodynamics shows that states having 
minimum entropy production compatible with the system constraints are sta- 
tionary states 11 . In reality, this property characterizes the case of linear flux 
laws with constant phenomenological coefficients. Nevertheless, MREP allows 
for many non-equilibrium physical processes to be considered within a general 
framework 11 . Given the importance of this matter, it is convenient to recall 
some essential aspects of the non-equilibrium theory, leaving details to the ded- 
icated treatises 12 . 

The basic equation of non-equilibrium thermodynamics is derived from the 
one by Gibbs. In general, the entropy change of a system can be written as 



where d e S is due to interaction with the system surroundings (actually, we are 
dealing with a closed system) and diS is the entropy production due to internal 
change of the system . In the case that internal changes are due only to chemical 
reactions , the entropy production term can be written as 11 



where /ik is the chemical potential of the k-component of a mixture contain- 
ing c chemical species and dn^ is the corresponding molar change. For the 
ideal system [ik = (k + RT In Nk where (k is a quantity independent of actual 
composition and Nk is the molar fraction of the k— component. If the sys- 
tem holds r chemical reactions, the molar change of the k-component can be 



A. Theory 



dS = d e S + diS 



(6) 




(7) 
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written as dn^ — Y^i=i ^kid^i where V)~i is the stoichiometric coefficient of the 
k-component in the /-reaction which shows a "displacement" d£i . Thus , the 
entropy production can also be written as 

diS = --= ^fikmdZi = = ^Tid^i (8) 

kl I 

where T/ = -J^k^kl = -J^k^klCk - R T J2k inN k M is tne chemical affin- 
ity of the /-reaction. At equilibrium the affinities vanish, so that vuCk — 
—RT'Y^j.YaNfc 1 . Thus, affinities can also be written as T; = — _RTlnIIfc (nt/rik) 
If at any instant the entropy changes as a function of chemical composition as 
well as other quantities characterizing the system, it is possible to write an 
equation for the rate of entropy production, that is, dS/dt = d e S/dt + diS/dt, 
where 

v i — d^i/dt standing for the flux or velocity of the /-reaction. It is assumed 
that near equilibrium fluxes are linear with respect to the affinities, that is, 
vi = Li m T m where L; TO are called the phenomenological coefficients 11 . 
As shown by Onsager 13 , based on the time reversal invariance of (microscopic) 
mechanical laws, the phenomenological coefficients form a symmetric matrix, 
that is Li m = L m i. It is to be pointed out that the criterion for sign assignments 
to stoichiometric coefficients is quite arbitrary. However, whatever the choice, 
if no external forces cause internal change, diS/dt ^ holds. 



B. The variational calculus 

The separation of entropy contributions given in eq. © may cause some 
problems of interpretation when the interaction with photons is considered. 
Classically speaking, the interaction with the e.m. field changes the density 
of internal energy with a rate given by the product E j where E stand for 
the strength of electric field and j for the vector of current density. Thus, 
a field releasing its energy to (emitted from) the system increases (decreases) 
the internal energy and thus the entropy. If configurational or chemical changes 
of the system are involved, the picture is slightly more complex. Actually, the 
e.m. field, which causes departure of the chemical reactions from equilibrium, 
reduces the entropy. A sequence of reactions brings the system to equilibrium, 
thus producing positive entropy. At the end of the process, the net change of 
chemical entropy is null. But reactions may produce heat and photons which 
alter the internal energy and thus the system entropy. Thus chemical reactions 
behave like a machine converting the absorbed photons to heat (and photons 
with lesser energies). About the entropy balance, further considerations will be 
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advanced in the Sec. IV. Now it appears convenient to write the rate of internal 
entropy production diS/dt as the sum of two contributions, that is, 

d^S , > 

—j£ = (?ext + first (10) 

where a ex t means the change of internal entropy due to external force and 
(Tint that due to the internal forces, that is, those bringing the system to 
equilibrium. The definitions of these two terms are to be searched for by means 
of the composition changes induced by absorbed photons, that is, 

r 

dn k = c k $dt + ^2 VkldZl (11) 
l=i 

where $ stands for the flux of the whole absorbed photons (in suitable units) 
and Ck for the fraction of photon flux inducing molar change of the k-component. 
It follows from eqs © and (|llf) that 



d 4 = -I*r. , 

dt T T 



where T$ = c^fJ-k m the following will be referred to as external affinity or 
"force" 11 . We define 

(Text = -y$r$ (Tint = ^ ^ I>/ . (13) 

I 

The equation (|12fl can be generalized to include more independent sets of chem- 
ical reactions by adding further index labelling reaction sets, that is, 



^E r ^- (14) 



Vint r£ 

ii 

Variational procedures allowing for the MREP must account for the con- 
straints that keep the system from equilibrium. These can be formalized by 
means of relations among affinities which are to be inserted in a Lagrange min- 
imization procedure, where the affinities are the functional variables. We will 
see in Sec. 3 that the constraints have the form 

r$-^r 7i = 7 = l,...,r. (15) 

3 

It is to be remarked that minimization concerns only a int , that is, only the 
internal processes bringing the system to equilibrium. However, to get a com- 
plete definition of the steady state, it is convenient to apply the minimization 
procedure to the function 

= 2<j ext + <ji n t + 4 At I r$ - r 7J j , (16) 
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by including T$ among the functional variables. Factor 2 takes into account that 
(Tint is a quadratic form on the affinities (this procedure is substantially different 
from that used in ref. [10]). As a result, it is obtained that $ = J2 7 ^7/2 and 
v 1 = v 1 i = X y /2 I = 1, r , that is, $ = ^ 7 u 7 . Thus, at the steady state the 
fluxes within a reaction set are the same for all reactions and the whole flux of 
reaction sets is equal to the flux of absorbed photons. 

In this paper we are concerned with reactions allowing for a diagonal On- 
sager's matrix (see also appendix A). The diagonal Onsager's matrix is peculiar 
to reactions that do not show interference effects so that 11 



v y i = LylTjl l = l,...,c. (17) 
Thus, by taking into account eq. I|15|l it follows 



(18) 
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Note that, at the steady state, flux is independent of the reaction index but 
depends only on the reaction-set index. The reaction affinity is related to the 
external force by the equation 

r 7i = r $ 1 -jf . (19) 

Now, by defining the RRs as 



R lX = 1/Lji , (20) 

we are able to describe a set of chemical reactions, activated by the photon flux 
with external force T$ , as a current u 7 passing through a series of resistances 
JL,j to which is applied a voltage T$. Thus, the drop in voltage on the resistance 
i? 7J - , that is, T 7 ; = i? 7 jU 7 , can be calculated by means of eq. (|19f) which gives 
the voltage partition, that is, 

r 7i = ^^-r $ . (21) 

Finally, for full correspondence to Ohm's law we must give a meaning to the 
Joule equation for the dissipated power, that is, 

W 1 = r$u 7 = u 7 2J r 7J - = T (Tint (22) 



which corresponds to the rate of entropy production of the 7-reaction set mul- 
tiplied by temperature. This has an evident meaning: as the electrical work 
dissipate the electric potential energy by producing heat, the chemical reactions 
dissipate the stored chemical energy by producing entropy (in general there is 
production of heat and photons) . 
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Variational calculus can be addressed to obtain information about the dis- 
tribution of reaction fluxes. To this end we should minimize Oi n t with respect 
to the fluxes by putting u 7 ; = u 7 I = 1,2, ...,c, T 7 ; = v^i/L^i — v^/L^i 
7 = 1, 2, r and $ — ^ 7 u 7 = 0. Minimization is to be applied to the function 



(23) 



that is, explicitly, 



where $ is to be included among the functional variables. By defining the 
equivalent phenomenological coefficient of the 7-reaction set as 

LyEQ = l/J2jT ( 2 5) 
I ll 

it is obtained from the variational calculus that r$ = x/2 and u 7 = L 7 bqx/2, 
that is, 

u 7 = L 1 eqT^, = <&L 1 eq/ E L 1 eq (26) 

7 

which is the equation of the partition current provided we define the equivalent 
resistance of the 7-reaction set as (see eqs. 1201 and 12*01) 

c 

R-lEQ = V L~ t EQ = 2J ^7J ' (27) 

III. APPLICATION TO ELECTRON-HOLE 
RECOMBINATIONS 



Until now, for simplicity we have considered cases where internal changes are 
to be ascribed to chemical reactions. To deal with electron-hole recombinations 
we should consider more suitable units , that is, densities (cm -3 ) rather than 
molar concentrations. To this end, we must divide the rate of entropy production 
by molar volume, Vm (cm 3 ), so that 

— = -77 fcTlnHfc (n k n k ) 

Vm at 

where Nq stands for the Avogadro number, the other symbols having the usual 
meanings. Now, cLNq^/Vm /dt is the reaction flux with the desired units, that is, 
cm _3 s -1 . In the following, as a consequence of this unit choice, the affinities 
will be calculated as 
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r = -fcTlnII fe (n k /n k ) 



A. The SK and LK cases 



Near equilibrium (but also in most practical cases) the rate equations can 
be written in the approximation of non-degenerate statistics, that is, emission 
rates of trapped carriers released into bands are independent of the occupation 
of band levels. In this approximation, rate equations assume a very simple 
form 8,10 . For simplicity's sake, it is convenient to begin by considering the cases 
of SK and LK . Formally, they can be dealt with as a single case. Indeed, 
they differ only in what kind of carrier is trapped and what is recombined, that 
is, if an electron is recombined with a trapped hole or, conversely, if a hole is 
recombined with a trapped electron. As stated in Section I, this difference has 
no formal relevance since it is due only to actual level position with respect to 
the equilibrium Fermi level 

Let N t mean the density of defects, p t the actual density of trapped holes 
(traps empty of electrons), n t the density of traps empty of holes ( filled with 
electrons), so that N t = p t + n t . Owing to interactions with ionizing radia- 
tions, electrons and holes are injected into conduction and valence bands respec- 
tively, with a rate Carriers are captured with probability c e , for conduction 
band electrons, and c p , for valence band holes, respectively. Thermal releasing 
of electrons into conduction band and of holes into valence band occurs with 
probabilities s e and s p (included is the interaction with black-body radiation), 
respectively. Thus, the net rate of electron capture is 8 

v n = c e np t - s e n t (28) 

and that of hole capture is 



v p = Spp t - c p pnt (29) 
At equilibrium v n = v p = 0. Thus, 

^ = * (30) 
n t c e 

m = s_p 

Pt Cp 

To define affinities, we need to fix a positive direction for processes as a sign 
reference for stoichiometric coefficients. Let us take as the positive direction 
that of the arrow pointing from the conduction to the valence band. Thus, the 
terms representing processes which bring electrons towards the valence band, as 
well as holes towards the conduction band, are associated with the stoichiometric 
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coefficient equal to +1. The opposite sign is associated with the terms describing- 
processes in the opposite direction. On this ground, the electron affinity is 



-fcTln 



(32) 



with v a = v c = +1, since n and pt appear in a term of positive direction, and 
Vb = —1, since n t appears in a term of negative direction. Thus, we can write 



fcTln 



np t n t 

The hole affinity can be obtained in an analogous way, that is, 



The external affinity is 



T„ = fcTln 



-fcTln 



Pt put 
pn t p t 



(33) 



(34) 



with c a = q, = — 1 since the reaction is opposite to the hxed positive direction. 
Thus 



r* = fcTln 



np 
Tip 



(35) 



so that the loop constraint l(TB)l is satisfied, that is, T$ — (I"i + T 2 ) = 0. For 
steady states near equilibrium we can be write, approximately, 



— , -, n t np t 

c e np t 1 — 

np t n t 



, P pnt 

v p = c p pnt 1 



pnt P t 

In the same approximation, the latter eqs. can be re-written as 



_ Cppnt 



irp — LpTp 



where 



L,. 



c e np t 



c p pn t 



kT p kT 

are the Onsager's coefficients. Finally, by definitions, the RRs are 

1 _ fcT 
L n c e Wp t 



(36) 
(37) 

(38) 
(39) 
(40) 
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L p c p pn t 

which agree with that found by SR (see eqs. 10 and (J3J). Note that the steady 
current is 



v = v n — v„ = 



iri n + Hp 

so that 

v (R n + R p ) = r$ . 

By taking into account the definitions of QFLs for free carriers (see eq|5Jl, it 
follows, 

r$ = kTlnnp/np = F n — F p , 

so that 

v (R n + R p ) — F n — F p , 

which is exactly the result found by SR (eq. EJ. 

Note also that , according to SR (by integrating the eq. 2.9 of ref. [8] over 
the whole states of conduction band) 

!p = ££ = Nn cxp [(E t - E C ) /kT] (41) 

where E t stands for the energy of trapping level. Now, by taking into account 
that 

- = -\ - exp(Et - F t ) (42) 
Pt 1 - ft 

ft being the occupancy of electron traps with QFL F t , it follows from eqs 
and that 

r„ = F n - Ft . (43) 
Analogously, for the affinity of trapped hole it can be shown that 

r p - F t - F p . (44) 

The eqs Q43JI and Q44J1 allow for a full correspondence to the SR results. 



B. The PW case 

SR in their paper considered only recombinations in single level centres, as 
in SK or LK processes. Now, it is advisable to consider also the case of two-level 
recombination centres as in the PW processes. With respect to the SK (or LK) 
two-step processes, in the PW case we must consider a further step, that is, the 
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recombinations of electrons trapped in levels labelled, say, 1 with holes trapped 
in levels labelled, say, 2. The corresponding recombination rate is 



= im lt p 2t - SuPunzt (45) 

where 7r stands for the recombination probability, n\ t and p\ t for the densities of 
electrons and holes trapped at the level labelled 1 respectively, n\ t andpit for the 
densities of electrons and holes trapped at the level labelled 2 respectively and 
s w for the probability of pair production by thermal excitation. At equilibrium 

ni^ = (46) 

P\t n 2t V 

so that 

{, Pitn 2t n lt p 2t \ 

= irn lt p 2t 1 • (47) 

V n lt p 2t Pu n 2t ) 

It is easy to verify that the affinity is 

r^fcTln* 1 *" 3 *^. (48) 

n\tP2t Pu n 2t 

Thus by using the definitions of electron and hole affinities (obtained in the 
previous section) suitably modified to account for the index levels , it is easy to 
prove that T$ = T„i + T w + T p2 . Now, the flux of the recombination process is 

= L^T^ = L w [T« - (r„i + I»] (49) 

where 

L u = Trn lt p 2t /kT (50) 

to which is associated the RR 

Ru, = 1/L U = kT/Trn lt p 2t . (51) 
The current through the recombination channel is 

v = ^ . (52) 

Note that the affinity T w is related to QFLs by the equation T w = F n — F t2 . 



C. The case of band-to-band recombination 

As a final example, let us calculate the RR associated with band-to-band 
recombinations which now appears as an easy task. Briefly, the net rate of 
recombinations is vq = iXQnp ~ sq where txq stands for the probability of 
electron-hole pair recombination and sq for the thermal emission of electrons 
from the valence to conduction band. At equilibrium rvp = sg/^g- Thus 
vg — nGnp(l — np/np) = LgTg where Tg — T$ and Lg — ncnp/kT. The 
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associated RR is Rg = ^/Lg = kT/iranp, so that the current can be written 
as v G = T$/R G . 



IV. FINAL REMARKS AND CONCLUSIONS 



At the steady state, the entropy of the system is time independent so that 
dS/dt = , that is, d e S/dt = —diS/dt . From eq. (|T2|l and by taking into 
account of the results of variational calculus we see that diS/dt = 0. This 
result is expected since the chemical composition is time independent as well. 
Nevertheless, it is worth dwelling upon this point to complete the considerations 
in Section II. Actually, the entropy balance is based on the following scheme: 

photons (his) — ► Chemical reactions — ► heat + photons (hv < hv) 

A flux of energy carried by photons, owing to internal reactions, is partly con- 
verted to heat and partly to photons (of lesser energy). The rate of heat 
production is thus dQ/dt — <f> (hv — hv') which, divided by the system tem- 
perature, is to be considered as the indirect contribution to d e S/dt due to the 
flux of ionizing photons. But the system is not isolated, so that the heat pro- 
duced by reactions is thus released to surroundings. This causes at the steady 
state d e S/dt = but a positive entropy is released to the surroundings thus 
contributing to the increase in the entropy of the universe. 

The definitions of RR rely on the diagonal property of the Onsager's ma- 
trix associated with the recombination processes. Besides some nice cases, the 
generalization to non-diagonal processes does not appear to be possible. In Ap- 
pendix B we try to diagonalize the thermodynamic problem of a two-reaction 
process that satisfies eq. <|2U|) . 

We again stress that the results found in this paper hold in the linear range of 
non-equilibrium thermodynamics. In this connection, it is understood that they 
are to be considered in the same spirit as one considers the linear approximation 
when small deviations from linearity can be disregarded . The validity of 
MREP far from equilibrium is still debated 15 . Moreover, in general, relations 
among fluxes and forces are not linear. However, linearization of relations is 
possible provided the system dealt with is near the steady state considered 
as the "reference state" 16 . Unfortunately, the matrix of coefficients is not 
symmetric for long. Thus, at the moment, definitions of RRs (which conserve 
thermodynamic meanings) for processes far from equilibrium appear a difficult 
task. 

In this paper the electron-hole recombinations in a excited material near 
equilibrium has been investigated within the framework of the non-equilibrium 
thermodynamics ruled by the MREP principle. A formal Ohm-like law has 
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been proved to hold for recombinations, thus leading in a natural way to the 
concept of recombination resistances, as defined by Shokley and Read. The 
entropy produced by recombinations is found to obey a Joule-like law where it 
(multiplied by temperature) play the role of the " dissipated power" . This makes 
the correspondence between steady electron-hole recombinations and electrical 
processes complete, at least when linear flux laws hold. 

APPENDIX A 



The reactions dealt with in this paper can be represented by means of a 
chain-process scheme, that is, 

n 1 ^±n 2 ^ ± ... <=* rii <=* ... ^ n c 

where the I — th reactions depend only on the I and the I + 1 components. This 
picture leads to a simple form for affinities, that is, 

with y^i+i = — vi+x,l = 1 • By definition it follows that vi = LiTi thus giving a 
diagonal Onsager's matrix. Now, let us suppose that the set of reactions form 
a closed loop, that is, 

n 2 <=► ... <=> n c _i 

TI Tl 

rii <=* n c 

In this case J^i — 0' that is 

^ T; = -T c = - (v Ct in c + fi, c /ii) = H\ - n c 

If we consider an external force that moves the system from equilibrium in the 
point between components 1 and c, then it is opposite to the r c force that leads 
the system towards equilibrium. Thus, by labelling this force with $ we can 
write 

r* = -r c = ^r ; 

l^c 

from which we obtain eq. I)15[l. 
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Now let us consider a case where eq. l(T5t holds but the symmetric Onsager's 
matrix is not diagonal. Let us consider a two-component reaction set with fluxes 

v i = LuTi + L12T2 
v 2 = L21T1 + £ 22 r 2 

which can be re-written as 



V\ = (in - L12) Tx + Li 2 r$ 

l>2 = L 2 lT$ + (L22 — L21) T2 



At the steady state v\ — V2 so that 



(ill —V) + i L 22 - V) 

where y£ i and r\ = L12 = L 2 \. Moreover we see that 

(in - 77) (L22 - rj) 



(in - 77) + {L22 - rj) 



1 



which does not has the formal structure of Ohm's law. But, by defining the new 
fluxes 



(in - 77) (L22 - rj) 



(in - 77) + (L22 - 77) 
and the new phenomenological coefficients 

La = La — 77 

we are able to diagonalize the thermodynamic problem provided a new flux is 
considered, that is, 

Consequently, the variational calculus is to be applied to the function 



n 



v\T x u 2 r 2 u$r$ 1 



a [r* - (ri + r 2 ) 



fjr\ rj~i rjn rj~i rjn 

Now, we can define the reaction resistances as Ri = 1/ (La — 77) and R$ = \ jr\ 
, where i?$ is to be considered in parallel with the other two resistances. Note 
that in the case = 77 it follows that V\ = V2 = 77r$, and v\ = w 2 = which 
is consistent with the infinite resistance. 
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